Ground-state phase diagram of the one-dimensional dimerized t—J model 

at quarter filling 
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The ground state of the one-dimensional dimerized t—J model at quarter filling is studied by a 
Lanczos exact-diagonalization technique on small clusters. We calculate the charge gap, spin gap, 
binding energy, Tomonaga-Luttinger-liquid parameter, Drude weight, anomalous flux quantization, 
etc. We thereby show that the two types of dimerization, i.e., a dimerization of hopping integral 
and a dimerization of exchange interaction play a mutually competing role in controlling the ground 
state of the system and this leads to the emergence of various phases including the Mott insulating, 
Tomonaga-Luttinger-liquid, and Luther-Emery-liquid phases. The ground-state phase diagram of 
the model is given on the parameter space of the dimerizations. 
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I. INTRODUCTION 

There are a number of low-dimensional correlated elec- 
tron systems described by the Hubbard and t—J mod- 
els at quarter-filling. One of the examples is the Bech- 
gaard salts (TMTSF)2X and (TMTTF)2X with X^PFg, 
CIO4, Br, etc., where there are three electrons in the 
two highest-occupied molecular-orbitals of a dimerized 
molecule, e.g., (TMTTF)2, and the system is at | filling 
in terms of electrons, which corresponds to the quarter 
filling in terms of holes In this system, dimerization 
of the molecules is known to play an essential role: de- 
pending primarily on the strength of dimerization, there 
appears a variety of electronic phases, such as antiferro- 
magnetic insulating, paramagnetic metallic, and super- 
conducting phases Some theoretical calculations of 
the dimerized Hubbard models have been done to clarify 
the nature of this system 

Another example is the transition- metal oxide NaV2 05 
which is reported to be a quarter-filled ladder system, 
exhibiting a spin-Peierls-like phase transition accompa- 
nied by a charge ordering jj]. This system may also be 
regarded as a dimerized system at quarter filling if we 
may assume that two V ions on each rung of the ladder 
form a dimer molecule Q . Effects of the lattice dimeriza- 
tion on the correlated electron systems in two dimension 
(2D) have also been studied in connection with cuprate 
superconductivity, where the spin-gap phenomena above 
Tc have attracted much attention. The simplest possible 
model that exhibits a spin gap which can survive against 
hole doping may be a dimerized t—J model with an al- 
ternating exchange interaction, where an enhancement 
of the singlet superconducting correlation has been sug- 
gested §. 

Motivated by such developments in the field, we study 
in this paper the one-dimensional (ID) dimerized t—J 
model at quarter filling, of which not very much is known 



so far. The model we study is defined by the Hamiltonian 

H=-hY, i4aCja + H.C.) - t2 ^ {cijla + H.C.) 
{ij)cr {kl)a 

.-^) (1) 



J2^(Sfe-S 



{kl) 



where 



cl^=c\^{\—ni-a) is the constrained electron- 
creation operator at site i and spin a {=],[), Si is the 
spin-i operator, and rii is the electron-number operator; 
hereafter we refer to the fermionic particle as 'electron', 
which corresponds to, e.g., the hole in the organic com- 
pounds. We use the ID lattice shown in Fig. 1; {ij) 
stands for nearest-neighbor bonds with parameters ti and 
Ji and {kl) for those with parameters ^2 (^^i) and J2 
(>Ji). The model tends to the usual t—J model when 
there is no dimerization (^1=^2 and Ji=J2), of which 
much study has been made |[^-[lo[| , whereas in the limit of 
strong dimerization, the model represents an assembly of 
isolated dimers. We retain the relations between param- 
eters t and J obtained from perturbation, i.e., Ji—At\/U 
and J2—At^/U , in order to reduce the number of param- 
eters, where U is the corresponding on-site Hubbard in- 
teraction. We thereby keep a relation Ji/ J2={ti/t2f ■ 
We thus have three independent parameters, and if we 
take t2 as a unit of energy, then we have two parame- 
ters, for which we will take parameters representing t- 
dimerization and J-dimerization (of which a specific def- 
inition is given below). 

We employ a Lanczos exact-diagonalization technique 
jm , which is used to obtain energies of the ground state 
and a few low- lying excited states. We denote the number 
of lattice sites by Ng and the numbers of up- and down- 
spin electrons by iV| and iVj^, respectively. The electron 
density is then given hy n— {N^ -f Ni)/Ns, so that n—1 
represents the half-filled system. Here we restrict our- 
selves to the case of n=l/2, i.e., quarter filling. We use 
the finite-size systems of the size 4, 6, and 8 unit cells (or 
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8, 12, and 16 sites, respectively). In order to achieve a 
systematic convergence to the thermodynamic limit, we 
choose periodic boundary condition for N^ + Ni = 4m +2 
and antiperiodic boundary condition for + = Am 
where m is an integer. 

We will examine the ground-state properties by calcu- 
lating the charge gap, spin gap, Luttinger liquid param- 
eters, anomalous flux quantization, etc., as a function of 
the dimerization strength. We will thereby show that 
the competition between the two types of dimerization 
leads to various ground-state phases such as a Mott insu- 
lating phase, Tomonaga-Luttinger-liquid phase, spin-gap 
phase, etc. The ground-state phase diagram of the model 
is thereby given. 

This paper is organized as follows; In Sec. 11, we 
present the calculated charge and spin gaps and clar- 
ify the mechanism of insulator-metal (or superconduc- 
tor) transition. In Sec. Ill, we calculate the Tomonaga- 
Luttinger-liquid properties of the model and also discuss 
a possibility of singlet superconductivity. In Sec. IV we 
present a phase diagram of the system by summarizing 
the results given in Sec. Ill and IV. Conclusions are given 
in Sec. V. 



II. INSULATOR-METAL TRANSITION 

In this section, we calculate the charge gap, spin gap, 
and binding energy of the model, and discuss the mech- 
anism of insulator-metal (or superconductor) transition 
of the system. First let us introduce two parameters for 
dimerization; i.e., the strength of dimerization in the hop- 
ping integral defined as 



t2_tl 
tl 



(2) 



which we call i-dimerization, and the strength of dimer- 
ization in the exchange interaction defined as 



Jd 



•h-Ji 
tl 



(3) 



which we call J-dimerization. We also take J2/i2 as a 
measure of the strength of J-dimerization because if we 
keep td constant (>0) then J2/t2 oc Jd- These are the 
key parameters which control the electronic state of the 
system; the i-dimerization has the effect leading to the re- 
pulsive interaction among electrons that acts when differ- 
ent spins come in a single dimer, and the J-dimerization 
has the effect promoting the spin-singlet formation be- 
tween spins coming in a single dimer. These effects man- 
ifest themselves in the electronic state of an isolated t—J 
dimer with an electron. The single-particle gap of the 
dimer is given as ?7dimcr=2i— J, where t and J are the 
hopping and exchange parameters of the single bond, and 
thus C/dimcr may be regarded as the effective Hubbard- 
like interaction of the dimer. The value of C/dimcr can 



be either positive or negative depending on whether the 
value J/t is smaller or larger than 2. The competing 
effects of the two types of dimerization may thus be ex- 
plained in the limit of strong dimerization. Now let us 
examine whether the competition between the effects of 
two types of dimerization in the ID system can lead to 
the insulator-metal (or superconductor) transition. For 
this purpose we calculate the charge gap, spin gap, and 
binding energy. 



A. Charge gap 

The charge gap may be defined by 
2 



Eo{N^ + \,N{)-E^{N^,NC)\ 
-[Eo{N^,N^)-Eo{N^-l,N^)\ 



(4) 



where Eo{N^, N^) is the ground-state energy of the sys- 
tem with Nij up-spin and down-spin electrons. This 
is evaluated for a number of clusters; throughout the pa- 
per we use a linear extrapolation of the data obtained 
from finite-size clusters with respect to 1 /Ng to estimate 
the value at the infinite system size because to use more 
sophisticated scaling formula will give no qualitative dif- 
ference to the obtained results although unphysical small 
negative values sometimes appear due to the errors of this 
scaling. The calculated results for Ac as a function of the 
t- and J-dimerizations (i.e., td and J2/t2) are shown in 
Fig. 2 (a). 

In case where there is no dimerization (ii=t2 and 
Ji=J2), i.e., for the homogeneous t—J model, a num- 
ber of studies [^^^Oj have shown that there is no charge 
gap Ac/t2=0 in the entire region of J2/t2 (except in the 
region of phase separation). Our result is consistent with 
this. In case where there is dimerization, the charge gap 
opens as seen in Fig. 2 (a). At J2/t2~^Q the system is 
equivalent to the ID dimerized Hubbard model at U—>oo, 
so that the dimerization gap of the size Ajj = 2{t2—ti) 
opens at the Brillouin zone boundary, which is the charge 
gap. In other words, the separation between the low- 
est unoccupied (antibonding) band and the highest com- 
pletely occupied (bonding) band is given by . In real- 
space picture, there is one electron per dimer, i.e., the 
number of electrons is equal to the number of dimers, 
and by regarding the dimer as a site one may have an 
effective half-filled band with the effective Coulomb re- 
pulsion of [/off — Ao. We thus have a Mott insulator due 
to i-dimerization. 

With increasing J-dimcrization J2/t2, we find 
that Ac/t2 decreases and becomes zero at a value 
J2/i2=(J2/t2)c''"S'' where the gap closes. (J2/i2)f 
represents the critical strength of J-dimerization at which 
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a transition from insulating phase to metallic (or su- 
perconducting) phase occurs. This may be explained 
as follows; In the insulating phase (Ac/i2>0), the ef- 
fective repulsion t/gff is given by the difference between 
loss of the kinetic energy and gain in the exchange in- 
teraction when one brings two electrons into a dimer, 
i.e., UeS='^{t2—ti)—J2- Thus, with increasing J2/i2, the 
measure of the gap UeS decreases, and at some value of 
J2/i2 the charge gap closes. In the small J-dimerization 
limit, we may use the perturbation theory with respect to 
J2/t2 (or Ji/ti) 1^. After a straightforward calculation 
we obtain the result for the charge gap: 



2(t2-ii) 



^ In2/Ji ^ J2\^^ 4(^i+t2) 

TT \tl 12) t2 -tl 



(5) 



In the limit of J2^0 and Ji— >0, this reduces to Ac=2(<2 — 
ti), which is equal to defined above. We compare 
the exact-diagonalization data with this analytical ex- 
pression in Fig. 2 (b) where we find a good agreement. 

With increasing J2/t2 further, we find that at some 
J2/i2 value the charge gap starts to increase again (see 
Fig. 2 (a)), at which the spin gap opens as we will see 
below. This is because Ac reflects the effect of singlet 
binding energy in the charge-gapless region. 



B. Spin gap and binding energy 

The spin gap may be defined by 

A,^Eo{N^ + l,N^^l)-EoiN^,N^l). (6) 

The calculated results for As as a function of the t- 
and J-dimerizations are shown in Fig. 3. When there 
is no dimerization, the results obtained are consistent 
with a recent result for the t—J model [|l^. In a small 
J2/t2 limit we have A^— >0 because there is no spin cor- 
relations in the system. With increasing 72/^2 with a 
fixed strength of t-dimerization, we find that the spin gap 
opens at some J2/i2 value, which we define as (^2/^2)0^'", 
the critical strength of J-dimerization at which the spin 
gap opens. We find that the spin gap remains finite in the 
region between J2/t2 = {J2 /t2)f^^'^^^ and the larger J2/t2 
value at which the phase separation occurs. We also find 
that with increasing td, the critical strength (J2/i2)c^'" 
becomes smaller and at the same time As/t2 increases. 
This means that the spin gap is enhanced by increasing 
the i-dimcrization. 

We note that the relation (72/^2)^" > (^2/^2)^''"^° 
always holds. This suggests that there exist two types of 
metallic region in the model; one is the phase where there 
are both gapless spin and gapless charge modes, and the 
other is the phase where there is a gap only in the spin 
excitation, which are the so-called Tomonaga-Luttinger 
(TL) region and Luther-Emery (LE) region, respectively. 
We will discuss this in Sec. III. 



A simple picture may be given to the case of strong 
J-dimerization where the spin-gap formation is ensured. 
When even number of electrons exist, in the lowest order 
perturbation of 12/ J2 two electrons with opposite spins 
always make a pair on the dimerized bond and gain the 
singlet formation energy 3J2/4. The electrons hop only 
as a pair tunneling through the virtual pair breaking in 
the forth order of ii/ Ji or t2/J2- The effective Hamilto- 
nian is of the form 



Hes = -i{slsi+2 + H.c), 



where 



t J 



(7) 



(8) 



' V2 

is the singlet Bose operator at site i (even number) and 



ti t 



1 "-2 



(9) 



is the effective hopping parameter of the boson. The 
Hamiltonian Eq. (7) may be derived in the same way 
as the effective Hamiltonian of the attractive Hubbard 
model in the strong-coupling limit HJl^ . 

We also calculate the binding energy defined as 

As = [Eo{N^ + 1, iVx - 1) - So(7Vt, TVx)] 

-2[So(7V^ + l,7Vi)-So(iVT,iVx)], (10) 

which is negative if two electrons minimize their energy 
by producing a bound state, and indicates a possible su- 
perconductivity. This value is meaningless unless the sys- 
tem is metallic or J2/i2>( J2/i2)?"*5°. In Fig. 4, we show 
calculated results for As as a function of J2/t2, where we 
find AsCiO for all values of td>0 until J2 / t2={J2 1 t2)T'^ 
is reached, and at this critical point A b starts to decrease 
suddenly. This mean that a bound state of two holes is 
formed in the entire region of the spin-gap phase and 
the binding energy is enhanced as J2/t2 is increased. We 
also note that at constant J2/i2 the value of Ab increases 
with increasing of td] i.e., the singlet binding energy is 
enhanced by the f-dimerization. 



III. LUTTINGER LIQUID BEHAVIOR 

It is well known that the ID interacting fermion sys- 
tems may be related to the Fermi-gas model in the 
continuum limit, where there are two different regimes, 
the Tomonaga-Luttinger (TL) regime and Luther-Emery 
(LE) regime. The essential difference between the two lies 
in the spin degrees of freedom; the TL region is charac- 
terized by the liquid phase with both a gapless spin mode 
and a gapless charge mode, whereas in the LE region the 
charge degrees of freedom is described by the TL liquid 
but the spin degrees of freedom have a gap. According to 
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the TL liquid theory [H^, various combinations of the 
parameters Kp and K„ describe the critical exponents of 
correlation functions of the system. In the absence of 
magnetic field, the dimerized t—J model is isotropic in 
spin space, so that K„~l holds, and we are left with the 
only parameter Kp. 

In the TL region, the spin and charge correlation func- 
tions show a power-law dependence as 



(11) 



D = 2VrK„ 



(17) 



and 



Equations (14)-(17) provide us with independent condi- 
tions on Kp, Vc, and D, which can be used to evaluate 
the TL-liquid parameter and to check the consistency in 
the TL-liquid relations. 



A. Tomonaga-Luttinger-liquid parameter 



i+rl 



2ik-pr j^Kp+K, 



r 



(12) 



respectively. We see that, for Kp<l, 2fcF-SDW or 2kp- 
CDW are enhanced and diverged, whereas for Kp>l, 
pairing fluctuations dominate. In the LE region, on the 
other hand, the spin gap opens and the 2A;f-SDW cor- 
relation decays exponentially. Because the contribution 
from spin excitations vanishes, the critical exponent of 
2fcF-CDW also changes and the asymptotic form is given 
by 



2ikt 



i+rl 



7' 



(13) 



Thus, in the region with Kp>l, the singlet spin correla- 
tion dominates over the 2fcF-CDW correlation. 

The relations between the correlation exponent Kp and 
the low-energy behavior of the model given below are use- 
ful for the evaluation of Kp-, First, the parameter Kp is 
obtained from the charge compressibility k and charge 
velocity Vc as 



Kp = -^VcK. 



The charge velocity may be determined by 



TV, 



— (-Ei,s=o — Eqi) 



(14) 



(15) 



where i?i,s=o is the energy of the lowest charge mode 
(measured from the ground-state energy Eq) at a neigh- 
boring k point. The inverse compressibility is given by 



1 _ 1 d'^Ep 
ii?K, Ns dv? 
2 

~ iV 



Eo{N^+2,N{)-E^[N^,N^)] 



[Eo[N^,N^)-Eo{N^^2,N^)\ 



(16) 



The parameter Kp is also related to the Drude weight 
D, the weight of the zero- frequency peak in the optical 
conductivity a^j , and may be obtained by considering the 
curvature of the ground-state energy-level as a function 
of the threaded flux Ill-|ll 



The calculated results for the TL-liquid parameter Kp 
are shown in Fig. 5 (a) as a function of J2/t2 for various 
strength of i-dimerization where the 8-, 12-, 16-site clus- 
ters are used although the cluster-size dependence of Kp 
is small. 

A limiting case of J2/i2"*'0 may be considered first. 
The ground state can be obtained by using the first-order 
degenerate perturbation theory around Ji=J2=0, where 
the wave function is the same as that in the U/t—^oo 
dimerized Hubbard model, i.e., a product of the state of 
spinless fermions describing the charge degrees of free- 
dom localized at each dimer and the state describing the 
spin system. The 2fcF-SDW correlation is dominant here: 
there is no chance for superconductivity. Note that even 
if the additional three-site terms are present, which are 
obtained through the Schrieffer- Wolff transformation to 
derive the dimerized t—J model from the dimerized Hub- 
bard model, they do not bring any change in the wave 
function in the first-order perturbation. Consequently, 
the dimerized t—J model in the J2/t2—^0 limit for any td 
values gives us the value of Kp^l/2. We however note 
that, at J2/t2=0 with a finite i-dimerization, the charge 
gap is positive Ac/t2>0 and the system is a Mott insu- 
lator. 

With increasing J2/t2, at any i^, Kp increases, and 
in the intermediate strength of J2/t2 (and for a rather 
small td), there appears the region with Kp>\ where the 
superconducting correlation is the most dominant. In 
the g-ology, this means the region of 52 <0 is effectively 
realized due to the attraction interaction caused by the 
J-dimerization. 

In Fig. 5 (b), we show the contour map for Kp on the 
parameter space {ti/t2, ^2/^2) where the contour lines are 
drawn by using a spline interpolation. We find that, for 
a fixed J2/t2, Kp decreases as t-dimerization increases. 
But as shown above the spin gap are enhanced by the 
increase of t-dimerization, so that the phase with both 
A<,>0 and Kp>l is realized in a reasonably wide range 
of the parameter values. We thus confirm that the sin- 
glet superconducting phase indeed exists. We also find 
that the calculated spin and charge correlation functions 
fit very well with Eqs. (11) and (12) as seen in Fig. 5 (b). 
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B. Drude weight and flux quantization 

The calculated results for the Drude weight defined 
in Eq. (17) are shown in Fig. 6. We find that, as t- 
dimerization increase, the dependence of D on J2/t2 be- 
comes stronger, and in the strong t-dimerization limit, D 
approaches the value at around the phase-separation 
point. The numerical technique used here is to thread the 
cluster ring with a flux $ and study the functional form 
of the ground-state energy with respect to the threaded 
flux, Eo{^). In general, Eo{^) consists of a series of 
parabola, corresponding to the curves of the individual 
many-body states i?„(<i>). This envelope exhibits a pe- 
riodicity of one in units of the flux quantum ^o=hc/e. 
The function Eq{^) (or the Drude weight and superfluid 
density) also yields information on the phenomenon of 
anomalous flux quantization; one may simply include the 
effect of a constant vector potential along the ring by the 
gauge transformation 



At p-ij<P 



where 



27r $ 



(18) 



(19) 



The existence of minima at intervals of half a flux quan- 
tum, which is the anomalous flux quantization, clearly 
indicates the existence of pairing. In Fig. 7, we show the 
calculated results for the flux quantization. We see that 
for td=0 (i.e., the uniform t—J model), the anomalous 
flux quantization is not observed. It occurs in the region 
where the binding energy is negative. At any strength of 
t-dimerization we find that the anomalous flux quantiza- 
tion occurs for appropriate strength of J-dimerization. 

In order to check the consistency of the TL-liquid re- 
lations, we compare the charge velocity obtained by two 
independent methods: one is from Eq. (15), and the other 
is from the relation 



(20) 



derived from Eqs. (14) and (17). The result is shown in 
Fig. 8, where we find that the reasonable consistency is 
indeed achieved. 



which occurs at around 72/^2—3 when td^O and extends 
to the large id region. The critical strength of ^2/^2 is 
almost independent to td- The phase separated region is 
determined by the condition k<0. When t-dimerization 
is dominant over J-dimerization, we find that the sys- 
tem becomes a Mott insulator, where the phase bound- 
ary is determined from the calculated values of the charge 
gap. On the other hand, when J-dimerization is domi- 
nant over t-dimerization, the system becomes metallic or 
superconducting. We note that there always appears the 
region of the TL-liquid phase between the region of the 
Mott insulating phase and the region of the LE-liquid 
phase. The phase boundary between the TL-liquid and 
LE-liquid regions is determined from the calculated val- 
ues of the spin gap although whether the spin-gap region 
exists in the homogeneous t—J model at quarter filling 
is not clear from our exact-diagonalization data ||ic| ]. In 
the limit of strong t-dimerization (i.e., td^oo), the model 
represents an assembly of isolated dimers and the system 
is an insulator in the entire region of J2/t2 unless the 
phase separation occurs. 

As for a possible correspondence with experiment, we 
may refer to a Bechgaard salt (TMTTF)2X where the 
dimerization strength of td— 0-1 is reported and the Mott- 
insulator to metal transition induced by pressure has 
been observed |2^j2^. We could argue that our phase 
diagram includes this phase transition at around td—O-i 
and J2/t2— 0.3 provided that the organic system can be 
described by the dimerized t—J model with a reasonable 
range of the parameter values J/t. 

We have also examined the phase diagram of the two- 
dimensional (2D) dimerized t—J model at quarter filling 
p2| and found that the phase boundary between the Mott 
insulating phase and the liquid phase has quite similar 
parameter dependence. This suggests that, irrespective 
of the spatial dimensions, the same mechanism discussed 
in Sec. II A operates in the present insulator-metal (or 
insulator-superconductor) transition in the quarter-filled 
dimerized t—J model. The main difference between the 
ID and 2D systems is that in ID the transition from the 
Mott insulator is to the TL-liquid region whereas in 2D 
it is to the singlet superconducting region. 



V. CONCLUSION 



IV. PHASE DIAGRAM 

By summarizing the results for a number of quantities 
obtained in the previous sections, we now draw the phase 
diagram of the ID dimerized t—J model at quarter-filling 
on the parameter space of t- and J-dimerizations. The 
result is shown in Fig. 9. 

When the value of J2/t2 (or the strength of J- 
dimerization) is very large, the system is phase separated, 



We have studied the ground state of the one- 
dimensional dimerized t—J model at quarter filling by 
using an exact-diagonalization technique on small clus- 
ters to calculate the ground-state properties such as 
the charge gap, spin gap, binding energy, Tomonaga- 
Luttinger- liquid parameter, Drude weight, anomalous 
flux quantization, etc. Thereby we have shown that the 
two types of dimerization, i.e., a dimerization of hop- 
ping integral (called t-dimerization) and a dimerization 
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of exchange interaction (called J-dimerization), plays an 
essential and mutually competing role in controlling the 
electronic ground state of the system; The t-dimerization 
has the effect leading to the repulsive interaction among 
electrons in a dimer and the J-dimerization has the ef- 
fect promoting the spin-singlet formation in a dimer. 
The resulting noticeable features we have obtained are 
the following: (i) The competition between t- and J- 
dimerizations induce the metal-insulator transition, (ii) 
There always appears the region of the TL-liquid phase 
between the region of Mott insulating phase and the re- 
gion of the LE-liquid phase, (iii) The spin gap and singlet 
binding energy are enhanced by the increase of dimeriza- 
tions. We have summarized the calculated results as the 
phase diagram in the parameter space of dimerizations. 

Finally we would emphasize that the correlated elec- 
tron systems at (and around) quarter filling indeed ex- 
hibit interesting properties as we have seen and further 
studies should be pursued in other quarter-filled systems 
from both theoretical and experimental sides. 
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FIG. 1. Schematic representation of the ID dimerized t—J 
model. The nearest-neighbor bonds have either ti and Ji 
(thin solid line) or t2 and J2 (bold line). The unit cells, each 
of which contains two sites, are indicated by dashed lines. 

FIG. 2. (a) Charge gap Ac/t2 as a function of J2/t2. (b) 
Perturbation estimate of the charge gap compared with ex- 
act-diagonalization data. 



FIG. 3. Spin gap As/t2 as a function of J2/t2- 



FIG. 4. Binding energy As/f2 as a function of J2/t2. 

FIG. 5. (a) TL-liquid parameter Kp as a function of J2/t2- 
(b) Contour map for Kp on the (ti/t2, ^2/^2) plane where 
dotted line is the boundary between the finite spin-gap and 
gapless regions (left panel) and distance r dependence of the 
charge and spin correlation functions where dotted lines are 
from Eqs. (11) and (12) (right panel). 



FIG. 6. Drude weight D as a function of J2/t2- 

FIG. 7. Energy difference Eo{(l>)—Eo{0) as a function of an 
external flux (j>. The t-dimerization td increases from the top 
panel to the bottom. 



FIG. 8. Charge velocity Vc as a function of J2/t2- The 
solid lines with symbols represent the value estimated from 
Eq. (15) and the dashed lines with symbols represent the value 
estimated from Eq. (20). 
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FIG. 9. Phase diagram of the ID dimerized t—J model 
at quarter flUing on the parameter space of t- and 

J-dimerizations. Dotted Unc separates the LE-hquid region 
from the TL-liquid region. Contour of constant Jd is shown 
by dotted Unes. 



